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Building on mode-coupling-theory calculations, we report a novel scenario for multiple glass transitions in a 
purely repulsive spherical potential: the square-shoulder. The liquid-glass transition lines exhibit both melting 
by cooling and melting by compression as well as associated diffusion anomalies, similar to the ones observed 
in water. Differently from all previously investigated models, here for small shoulder widths a glass-glass line 
is found that is disconnected from the liquid phase. Upon increasing the shoulder width such a glass-glass 
line merges with the liquid-glass transition lines, featuring two distinct endpoint singularities that give rise to 
logarithmic decays in the dynamics. These findings can be explained analytically by the interplay of different 
repulsive length scales. 
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In the field of glassy slow dynamics, many experiments and 
simulations have been inspired in recent years by the mode- 
coupling theory of the glass transition (MCT) Jl]]. The theory 
deals with density autocorrelation functions § q (t) with wave 
vectors q, and predicts their long-time limits f q . While in 
the liquid state f q = 0, the glass state is defined by f q > 0. 
MCT was first applied to the hard-sphere system (HSS) where 
a liquid-glass transition was identified J2l, and confirmed by 
experiments . In addition to liquid-glass transitions, for cer- 
tain interactions MCT also predicts glass-glass transitions: In 
this case an existing first glass state with f£ transforms into a 
second distinct glass state with > discontinuously. Such 
glass-glass transitions were predicted for the square-well sys- 
tem (SWS) where the hard-core repulsion is supplemented by 
a short-ranged attraction flSEl- In the SWS, the first glass 
state is driven by repulsion like in the HSS and the second 
glass state is driven by attraction. The competition between 
these two mechanisms is responsible for the emergence of 
glass-glass transitions. Such a line of glass-glass transitions 
extends smoothly a line of liquid-glass transitions into the 
glass state and terminates in an endpoint singularity. Close to 
the endpoint singularity the dynamics is ruled by logarithmic 
relaxation [7]. The predicted logarithmic decays were iden- 
tified in computer simulations and establish the relevance of 
endpoint singularities for the description of glassy dynamics 
JH A second dynamical anomaly predicted for the SWS 
concerns a reentrant liquid-glass line that causes melting by 
cooling |0> H 0] . This prediction of MCT was confirmed by 
computer simulation (10] and by experiments in colloidal sus- 



pensions 11 1LI12H . 

In this work we replace the attractive length scale in the 
SWS by a second repulsive length scale 8 of the square- 
shoulder system (SSS) as shown in Fig. Q] The SSS can be 
considered the simplest potential with two competing inter- 
particle distances; it is applied to describe properties of metal- 
lic glasses like cerium or cesium \13], micellar iflill and gran- 
ular materials |f5], silica iflrjll and water lfl7ll . 
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FIG. 1: Square-shoulder potential with control parameters packing 
fraction (p = Tt.pt/ 3 /6, shoulder height T = u^/k^T, and shoulder 
width 8 for particles of diameter d at density p. 

In the following, the glass-transition diagrams are calcu- 
lated from the singularities of the MCT functional QJJ] 



f q [vjk)= E v qJcp f k f p = -^- 



(la) 



k+p=ij 



with vertex V given by wave vector moduli q, k, p, the static 
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structure factor S q , and the direct correlation function c g : 

Vq,kpfkfp,yq,kp= z PSqSkS p {q-\kc k +pC p ]} 2 /q 4 . (lb) 

k+p=q 

Quantities S q and c q are given by the interaction potential. The 
discretization of the functionals is chosen like for the HSS and 
the SWS jHH with a number of wave vectors M = 600 and 
a wave-vector cutoff dq max = 80, the static structure factors 
of the SSS are obtained within the Rogers-Young (RY) ap- 
proximation | Jjl 2(1; further details of the calculations shall 



be found in a subsequent publication 12111 . For specific values 
of the control parameters packing fraction (p, shoulder height 
r, and shoulder width 8, Eq. ([T]i exhibits singularities where 
the f q change discontinuously indicating liquid-glass or glass- 
glass transitions. At these singularities one can define the so- 
called exponent parameter A, < 1 that fixes all other critical 
exponents of the theory (U. While typically at liquid-glass 
transitions, X rs 0.7, A, approaches unity at the endpoint of 
glass-glass transition lines signaling the emergence of loga- 
rithmic decay laws 0]. 

Figure|2|a) shows the glass-transition scenario for 8 = 0.13: 
The liquid-glass transition line (diamonds) starts at the HSS 
limiting value of q> c HSS = 0.5206 (dotted line for small F) for 
vanishing shoulder width F. For increasing shoulders and up 
to r rs 1.5, the curve exhibits a shift in the transition pack- 
ing fraction cp(r) to higher values - the glass initially melts 
upon cooling. This trend can be traced to the evolution of 
the static structure factor which becomes sharper and moves 
to lower wave vectors for higher F; the corresponding pair 
distribution functions show a higher probability for particles 
being at larger distances from each other. Hence, larger par- 
ticle separations weaken the cage and are compensated by 
higher densities at the glass transition, cf. 112 111 for more de- 
tails and a comparison to the melting-by-cooling phenomenon 
in the SWS. For shoulder heights from 2k^T to 3.5A:b7\ i.e. 
for 2 < r < 3.5, the glass-transition curve bends downwards 
and reaches the limiting value of the HSS for the outer core 
$hss = ( Phss/( 1 + § ) 3 = 0.3608 (dotted line for large F). 

If the packing fractions and repulsive strengths are in- 
creased beyond the liquid-glass transition line, one encounters 
an additional line of glass-glass-transition singularities (filled 
circles). Differently from glass-glass lines in all models in- 
vestigated previously, this additional line is located inside the 
glassy regime, disconnected from any liquid-glass transition 
line. It is bounded by two endpoint singularities (open cir- 
cles) where the additional discontinuity in the f q vanishes, A, 
reaches unity, and logarithmic decay laws emerge. 

When increasing the shoulder width 8 further, the glass- 
glass and liquid-glass transition lines move towards each other 
and start to merge for sufficiently high shoulders at around 
8 = 0.145. Figure|2|b) shows the situation for 8 = 0.15: From 
r = 2.0 to 3.5, the former glass-glass line now indicates a tran- 
sition from the liquid directly into the second glass state. For 
r < 2.0 and for F > 3.5, the formerly isolated glass-glass line 
crosses the liquid-glass line and extends into the glassy regime 




FIG. 2: Glass-transition diagrams for the SSS. Diamonds (♦) indi- 
cate liquid-glass transitions, filled circles (•) show glass-glass transi- 
tion points terminating in two endpoint singularities (o). (a) 8 = 0. 13. 
Dashed lines display the respective limits for hard spheres of di- 
ameter 1 and 1+8. The line of glass-glass transitions is discon- 
nected from the line of liquid-glass transitions, and moves towards 
and crosses it for larger shoulder width, (b) 8 = 0.15. Crosses (x) 
indicate the four states discussed in Fig. [3] relating the glass-glass 
transitions to features in the static structure. The dashed curve ex- 
hibits the liquid-glass transition line for a reduced wave-vector cut- 
off of dq max = 40 demonstrating that without contributions of higher 
wave vectors to the functional in Eq. Q} the lines of glass-glass tran- 
sitions are absent. 



as two glass-glass transition lines. On both ends, the end- 
point singularities retreat individually in a generic way which 
is similar to the case of the SWS |5j]. 

Figure |2jb) also exhibits the reentry phenomenon melting - 
by cooling for small F as discussed before for 8 = 0.13, but 
in addition a second reentry phenomenon is found between 
r = 3.0 and 3.5 where melting can also be induced by com- 
pression. The trends can again be traced back to the behav- 
ior of the static structure 12111 : Increasing the packing frac- 
tion from low values to around (p w 0.45, the pair distribution 
function shows an increased probability for particle contact at 
both inner and outer cores. However, for values larger than 
cp sw 0.45, the contact at the inner core grows at the expense 
of contact at the outer core; interparticle contacts at the outer 
core are suppressed by the high density and particles now col- 
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lide much more frequently at their inner cores. Since the cage 
represented by the inner core is too loose to trigger glassy ar- 
rest at such density, this is compensated by lower temperature 
(or equivalently by higher F) at the glass transition. 

For the SSS, the interplay between the hard and the soft 
core is manifest in changes at the (1) principal peaks signal- 
ing different interparticle distances, and (2) by a beating at 
larger wave vectors that is introduced by the oscillation fre- 
quencies q and (1 + 8)q, respectively. As explained above, 
changes in the principal peaks are sufficient to understand the 
overall behavior of the liquid-glass transition lines, including 
both reentry phenomena. In the following we demonstrate that 
mechanism (2) is responsible for the glass-glass transitions. 

It is known from the theory of Fourier transformations that 
the discontinuities of the direct correlation function in space, 
c(r), determine the large-wave-vector behavior of the direct 
correlation function c q which enters the MCT functional in 
Eq. ( TTbl I22I1 . For the square-shoulder potential the dominant 
singularities are inner and outer core which in leading order 
may be modeled by two independent hard-sphere diameters. 
To obtain a qualitative picture, regardless of the specific form 
of the closure relation, it is enough to consider HSS in Percus- 
Yevick (PY) approximation for which c q is known analytically 
1 23J ■ For large enough wave vectors its asymptotic behavior 



Cq , reads 
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Cj (q) for the inner core and c 2 (q) for the outer core are 
then given by c\ sy (q) =B\ cos(q)/q 2 and c%(q) = B2Cos(g[l + 
8\)/q 2 , where factors of (1 +5) in q 2 were absorbed into B 2 - 
This analytical PY result for the asymptotic decay describes 
the numerical solution of the RY structure factor very accu- 
rately as seen in the upper panel of Fig. [3] While for the re- 
spective limits of small and large hard spheres, c^ y (q) and 
c^ y (q), are relevant only individually, for parameters in be- 
tween, one expects an additive effect from both cores. As- 
suming the PY result for both cores we get for the combina- 
tion c* s y(q)=cf y (q)+C2 y (q): 



c{q) 



asy 



= ^{2fi 2 cos([l+8/2]^)cos( 9 8/2) 

q +(£?!- fi 2 )cos(<7)} 
= \ {2Bi cos ([1 + 5/2]^)cos ( 9 8/2) 

" +(B 2 -Bi) cos ([! + %)} 



(3) 



For B\ = B2, the second lines in Eq. (0) vanish and a simple 
beating is obtained, for B\ ^ B2 additional terms remain. Fig- 
ure[3]demonstrates the evolution of the tail in c q for increasing 
shoulder height. Starting from the HSS of the inner core (up- 
per panel of Fig. [3]), we get the original 1/q 2 decay as inferred 
from Eq. (fJI or equivalently the second line in Eq. 0. Switch- 
ing on the shoulder, B\ > Bi > 0, the original tail is lowered 
to B\ — Z?2 and an additional modulation of Aq is obtained for 
the term with mixed frequency. For the amplitude variation 
of the beating one obtains dAq = 2k/ 8, so for 8 = 0.15 one 




FIG. 3: Direct correlation function q 2 c q for 8 = 0.15, (p = 0.50, and 
T = 0.0, 1.0, and 4.0, respectively from top, T = 4.0 and 9 = 0.35 in 
the bottom panel, cf. crosses in Fig.|2jb). Full horizontal lines dis- 
play the value of B((p) in Eq. l(2j for (p = 0.50. The dashed curve in 
the upper panel shows the asymptotic solution of Eq. lO, the dashed 
curve in the bottom panel shows the asymptotic solution of Eq. ifJJ 
with q replaced by (1 + 8)151. The emergence and vanishing of a beat- 
ing from top to bottom explains the occurrence of the glass-glass 
transition lines in Figs.[2]a) and|2lb) through additional contributions 
to the functional in Eq. (fTJ. 



gets Aq ~ 42d~ 1 , a beating of amplitude variation of qd w 42 
as seen in Fig.[3]in the two middle panels. If both diameters 
are of equal importance, the beating terms in the first lines of 
Eq. (0 become dominant as seen in the third panel of Fig. [3] 
If the outer core is dominant, the second equation in Eq. (0) 
is relevant: Now the decay for the larger core is supplemented 
by the same beating but with prefactor B \ . For yet higher con- 
tributions from the outer core the beating vanishes again as is 
seen in the bottom panel of Fig. [3] 

Figure[3]explains the origin of the glass-glass transition line 
by the evolution of additional contributions to the MCT vertex 
in Eq. ( fTbl ). While the upper panel shows a tail in c q that 
is small enough to no longer influence the MCT vertex, the 
subsequent panels show contributions to the integrals from the 
beating at wave vectors from qd = 20 to 80 that cannot be 
neglected. To show that it is these contributions that cause 
the glass-glass transition line we eliminate these contributions 
by shifting the wave-vector cutoff from qd = 80 to qd = 40 
and perform additional calculations. The result is shown in 
Fig. |2jb) as the dashed line. With the beating switched off, 
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the dashed line of liquid-glass transitions is recovered without 
any indication of crossings or glass-glass transitions. Figure[3] 
also demonstrates that the beating is only relevant in a finite 
region that is bound from above and from below which relates 
to both endpoints of the glass-glass transition line. 

To make contact with experimental systems, we redraw the 
transition diagram of Fig. |2jb) in a pressure-vs-temperature, 
P-T, diagram in Fig.@] using the RY thermodynamically con- 
sistent equation of state. For a path of constant T and vari- 
able P, the diffusivity of the dynamics varies with the distance 
from the liquid-glass transition line. E. g., for T w 0.35mo/£b 
and starting from low P, the diffusivity first decreases un- 
til P w 15Mo/flf 3 , then it increases anomalously until around 
P w 30mo/^ 3 , and then it decreases again. Such behavior, 
known as diffusion anomaly, is experimentally observed in 
liquid water l24ll . 



trol parameter, either T, (p or P). The relevance of the glass- 
glass transition scenario found here goes beyond the specific 
SSS potential: The competition between two repulsive length 
scales and the presence of diffusion anomalies are observed 
in several different contexts, ranging from soft matte r sy stems 
HSl to silica HH, and complex liquids JH E2I Si . ° ur 
work opens a perspective for understanding slow dynamics in 
these disparate systems. 
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Sp714/3-l and BMWi 50WM0741. EZ and FS acknowl- 
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and HES thank NSF Chemistry Division for support. 
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FIG. 4: Data from Fig.[2{b) for 6 = 0.15 shown with variables pres- 
sure P and temperature T. Crosses ( x ) delimit the approximate re- 
gion for a diffusion anomaly. Above and below that region, the dif- 
fusion coefficient decreases with pressure - within the region the dif- 
fusion coefficient increases with pressure P. 



In conclusion, MCT calculations predict that the square- 
shoulder system is characterized by an unusual slow dynam- 
ics, based on the competition between two repulsive length 
scales. This manifests itself in an unexpected glass-glass tran- 
sition line which can be completely disconnected from the liq- 
uid phase. Upon reducing the distance between the two repul- 
sive lengths, the glass-glass transition line moves towards the 
glass-liquid line, merging for a special value of the shoulder 
width and giving rise to logarithmic behavior of the decay of 
correlations in a wide region of temperatures and densities. 
The mechanism for the glass-glass transitions is derived ana- 
lytically from features of the static structure factor: The two 
length scales are reflected in a beating at large wave vectors in 
S q that cause additional contributions to the MCT functional 
and thus trigger an extra transition. 

The peculiar shape of the arrest lines echoes also in the liq- 
uid phase, generating diffusion anomalies (i.e. non mono- 
tonic dependence of the diffusion coefficient with the con- 
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